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I.   Introduction 

This  report  contains  a  target  position  estimation  program  for 
each  of  the  following  calculators:   The  Hewlett-Packard  HP-41CV, 
the  Sharp  PC-1500  (or  Radio  Shack  TRS-80  PC-2) ,  the  Sharp  PC-1211 
(or  Radio  Shack  TRS-80  PC-1) ,  the  Casio  FX-702P  and  the  Texas 
Instruments  TI-59. 

The  programs  provide  a  means  of  implementing  a  position  esti- 
mation procedure  that  is  described  in  Appendix  1.   The  procedure 
is  based  on  the  following  assumptions:   Bearings  taken  on  or  from 
a  target  are  available  from  two  or  more  stations  of  known  position 
The  positions  of  the  stations  and  the  target  are  such  that  they 
can  be  considered  to  be  on  the  surface  of  a  plane  (a  flat  earth 
model) .   The  error  in  a  bearing  taken  on  or  from  a  station  is  a 
normal  random  variable.   Its  standard  deviation  (bearing  error) 
is  known  and  its  mean  (bias)  is  zero  (if  there  is  bias  in  a  bear- 
ing, it  is  known  and  it  is  removed) .   The  errors  in  bearing  mea- 
surements are  independent. 

The  position  estimation  procedure  requires  an  initial  estimate 
of  the  target's  position.   In  the  programs  that  are  given  here, 
the  initial  estimate  is  at  the  intersection  of  the  bearing  lines 
determined  by  the  first  two  target  bearings  that  are  input  to 
the  program.   Because  of  the  use  of  this  method  to  choose  an  ini- 
tial estimate,  the  first  two  stations,  with  respect  to  the  order 
of  data  input,  should  be  chosen  so  that  the  intersection  of  their 
bearing  lines  is  likely  to  be  closer  to  the  target's  position 
than  the  intersection  of  the  bearing  lines  from  any  other  pair 
of  stations. 


The  target  ranges  and  errors  of  the  bearings  of  the  first  two 
stations  determine  the  distance  between  the  intersection  of  their 
bearing  lines  and  the  target's  position.   In  particular,  if  the 
angular  separation  between  two  stations  as  seen  from  the  target 
is  small  relative  to  the  bearing  error  of  one  or  both  of  the  sta- 
tions, the  bearing  lines  from  the  two  stations  may  not  intersect. 
If  they  do  not  intersect  and  they  are  not  parallel,  their  reci- 
procal bearing  lines  will.   In  this  case,  if  the  observed  bearings 
from  the  two  stations  were  the  first  two  bearing  inputs  to  the 
program,  the  initial  estimate  of  the  target's  position  would  be 
at  the  intersection  of  the  two  reciprocal  bearing  lines,  and  a 
gross  error  in  the  final  position  estimate  could  result. 


II.   Program  User  Instructions 

Before  using  a  program  for  the  first  time,  refer  to  the  notes 
that  follow  the  user's  instructions.   Also,  note  the  comments  in 
Section  I  regarding  the  relationship  between  the  accuracy  of  the 
computed  estimates  and  the  order  of  data  entry. 

For  a  PC-1211,  use  the  PC-1500  User  Instructions.   Use  the 
DEF  mode  and  substitute  SHFT  for  DEF  wherever  DEF  appears  in  the 
instructions. 

Station  positions  are  determined  with  respect  to  a  known  ref- 
erence point.   The  reference  point  can  be  a  station  position,  in 
this  case  the  station  bearing  and  range  from  the  reference  point 
are  both  zero. 

Program  listings  are  given  in  Section  IV  and  two  example 
applications  are  discussed  in  Section  III. 


Position  Estimation  Program 
HP-41CV  User  Instructions 


Step 


Instruction 


Prompt 


Press 


10 


11 


12 


Select  USER  mode  and  enter  the  program 
(see  Note  1) .   Press  Z+  to  run  the 
program. 

Key  in,  in  decimal  degrees,  the  obser- 
ved bearing  of  the   target  from  a 
station  or  the  reciprocal  of  the  ob- 
served bearing  of  the  station  from 
the  target  (see  Note  2) . 

Key  in  the  station  bearing  in  decimal 

degrees  from  the  reference  point.   Use 

zero  if  the  station  is  the  reference 
point. 

Key  in  the  station  range  (in  any  units) 
from  the  reference  point.   Use  zero 
if  the  station  is  the  reference  point. 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

Repeat  Steps  2,  3,  4  &  5  for  one  or 
more  additional  stations. 

Compute  bearing  &  range  estimates. 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

A  computed  target  range  estimate  in 
meters.   For  an  elliptical  contain- 
ment region,  go  to  Step  10  or  Step  17. 
To  enter  additional  data  from  new  or 
old  stations,  go  to  Step  24. 

For  a  containment  ellipse  of  a  given 
containment  probability,  press  LN. 

Key  in  the  desired  containment  proba- 
bility. 

Computed  value  of  the  ellipse  size: 
(See  Note  3. ) 


OBS  BRG? 


STA  BRG? 


STA  RNG? 


BRG=est. 


RNG=est. 


PRB? 


SIZE=val 


E  + 


R/S 


R/S 


R/S 


R/S 


1/x 

R/S 


LN 

R/S 

R/S 


Step 

Instruction 

Prompt 

Press 

13 

Computed  semi-ma jor  axis  length: 

SMJ=val. 

R/S 

14 

Computed  major  axis  direction: 

DRC=val. 

R/S 

15 

Computed  semi-minor  axis  length: 

SMI=val. 

R/S 

16 

Computed  containment  ellipse  area: 
(See  Note  4. ) 

A  =  val. 

17 

For  a  containment  ellipse  of  a  given 
size,  press  LOG. 

LOG 

18 

Key  in  the  desired  containment  ellipse 
size  (see  Note  3) . 

SIZE? 

R/S 

19 

Computed  value  of  the  containment 
probability: 

PRB=val. 

R/S 

20 

Computed  semi-major  axis  length: 

SMJ=val. 

R/S 

21 

Computed  major  axis  direction: 

DRC=val. 

R/S 

22 

Computed  semi-minor  axis  length: 

SMI=val. 

R/S 

23 

Computed  containment  ellipse  area: 
(See  Note  4.) 

A  =  val. 

24 

To  enter  additional  data  from  new  or 
old  stations,  press  /x  and  then  repeat 
Steps  2,  3,  4  &  5 

/x 

Notes : 

1.  The  program  size  is  35.   The  key  assignments  are:   TPE  -*■  E+, 
CON  ->  /x,  EST  ->  1/x,  SIZ  ■+   LOG  and  PRB  -*  LN.   If  they  are  not 
present,  they  must  be  made  in  order  to  follow  the  user 
instructions.   For  an  alternative,  see  Appendix  2. 

2.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

3.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 
probability  of  containment,  the  minimum  area  containment  region 
is  an  ellipse  centered  on  the  position  estimate.   The  semi- 
major  axis  =  ka„T  and  semi-minor  axis  =  ka.,T  where  k  is  the 

J  MJ  MI 

size  of  the  ellipse  and  o..^   and  a..-,    are  the  standard  deviations 

c  MJ       MI 

(uncertainty  measure)  of  the  position  estimate  in  the  major 
axis  and  minor  axis  directions. 

4.  The  area  units  are  the  range  units  squared. 


Position  Estimation  Program 
PC-1500  User  Instructions 


;tep 

Instruction 

Prompt 

Press 

1 

Enter  the  program.   To  run  the  program, 
press  DEF,  A.  For  a  PC-1211,  see  Note  5. 

DEF  A 

2 

Key  in,  in  decimal  degrees,  the  observed 
bearing  of  the  target  from  a  station  or 
the  reciprocal  of  the  observed  bearing 
of  the  station  from  the  target  (see 
Note  1) . 

OBS  BRG? 

ENTER 

3 

Key  in  the  station  bearing  in  decimal 
degrees  from  the  reference  point.   Use 
zero  if  the  station  is  the  reference 
point. 

STA  BRG? 

ENTER 

4 

Key  in  the  station  range  (in  any  units) 
from  the  reference  point.   Use  zero  if 
the  station  is  the  reference  point. 

STA  RNG? 

ENTER 

5 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

BRG  ERR? 

ENTER 

6 

Repeat  Steps  2,  3,  4  &  5  for  one  or 
more  additional  stations. 

7 

Compute  bearing  &  range  estimates. 

DEF  Z 

8 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

BRG=est. 

ENTER 

9 

A  computed  target  range  estimate  is  dis- 
played.  For  an  elliptical  contain- 
ment region,  go  to  Step  10  or  Step  17. 
To  enter  additional  data  from  new  or 
old  stations,  go  to  Step  24. 

RNG=est. 

10 

For  a  containment  ellipse  of  a  given 
containment  probability,  press  DEF,  X. 

DEF  X 

11 

Key  in  the  desired  containment  proba- 
bility. 

PRB? 

ENTER 

12 

Computed  value  of  the  ellipse  size: 
See   Note  2. 

SIZE=val. 

ENTER 

Step 


Instruction 


Prompt 


Press 


Computed  semi-major  axis  length: 

Computed  major  axis  direction: 

Computed  semi-minor  axis  length: 

Computed  containment  ellipse  area: 
(See  Note  3.) 

For  a  containment  ellipse  of  a  given 
size,  press  DEF,  S. 

Key  in  the  desired  containment  ellipse 
size  (see  Note  2) . 

Computed  value  of  the  containment 
probability: 

Computed  semi-major  axis  length: 

Computed  major  axis  direction: 

Computed  semi-minor  axis  length: 

Computed  containment  ellipse  area: 
(See  Note  3.) 

To  enter  additional  data  from  new  or 
old  stations,  press  DEF,  C  and  then 
repeat  Steps  2,  3,  4  &  5. 


SMJ=val 
DRC=val 
SMI=val 
A  =  val 


SIZE? 

PRB=val 

SMJ=val 
DRC=val 
SMI = val 
A  =  val 


ENTER 
ENTER 
ENTER 


DEF  S 

ENTER 

ENTER 

ENTER 
ENTER 
ENTER 


DEF  C 


Notes : 

1.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

2.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 

probability  of  containment,  the  minimum  area  containment 

region  is  an  ellipse  centered  on  the  position  estimate.   The 

semi-major  axis  =  ka_-_  and  semi-minor  axis  =  ka.,T  where  k  is 
J  MJ  MI 

the  size  of  the  ellipse  and  a„T  and  a.,_  are  the  standard 

^        MJ       MI 

deviations  (uncertainty  measure)  of  the  position  estimate 
in  the  major  axis  and  minor  axis  directions. 

3.  The  area  units  are  the  range  units  squared. 

4.  For  a  definition  of  the  program  initiating  keys  and  their 
function,  press  DEF,  H.   The  display  will  show: 

TPE  =  A  EST  =  Z  SIZ  =  S  PRB  =  X.   Next  press  ENTER.   The 
display  will  show:   CON  =  C.   To  repeat  the  display,  press 
ENTER. 

5.  For  a  PC-1211,  use  the  DEF  mode  and  substitute  SHFT  for 
DEF  wherever  DEF  appears  in  the  instructions. 


Position  Estimation  Program 

FX-702P  User  Instructions 

Step 

Instruction 

Prompt 

Press 

1 

Enter  the  program  (see  Note  1) .   To  run 
the  program,  first  press  Fl,  0. 

Fl  0 

2 

Next,  key  in  1  and  press  EXE. 

OPTION? 

1  EXE 

3 

Key  in,  in  decimal  degress,  the  obser- 
ved bearing  of  the  target  from  a 
station  or  the  reciprocal  of  the  ob- 
served bearing  of  the  station  from 
the  target  (see  Note  2) . 

OBS  BRG? 

EXE 

4 

Key  in  the  station  bearing  in  decimal 
degrees  from  the  reference  point.   Use 
zero  if  the  station  is  the  reference 
point. 

STA  BRG? 

EXE 

5 

Key  in  the  station  range  (in  any  units) 
from  the  reference  point.   Use  zero 
if  the  station  is  the  reference  point. 

STA  RNG? 

EXE 

6 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

BRG  ERR? 

EXE 

7 

Repeat  Steps  3,  4,  5  &  6  for  one  or 
more  additional  stations. 

8 

To  compute  bearing  &  range  estimates, 
first  press  Fl,  0. 

Fl  0 

9 

Next,  key  in  2  and  press  EXE. 

OPTION? 

2    EXE 

10 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

BRG=est. 

CONT 

11 

A  computed  target  range  estimate  is  dis- 
played.  For  an  elliptical  contain- 
ment region,  go  to  Step  12  or  Step  20. 
To  enter  additional  data  from  new  or 
old  stations,  go  to  Step  28. 

RNG=est. 

12 

For  a  containment  ellipse  of  a  given 
containment  probability,  first  press 
Fl,  0. 

Fl  0 

13 

Next,  key  in  3  and  press  EXE. 

OPTION? 

3  EXE 

10 

Step 

Instruction 

Prompt 

Press 

14 

Key  in  the  desired  probability. 

PRB? 

EXE 

15 

Computed  value  of  the  ellipse  size: 
See  Note  3 . 

SIZE=val. 

CONT 

16 

Computed  semi-major  axis  length: 

SMJ=val. 

CONT 

17 

Computed  major  axis  direction: 

DRC-val. 

CONT 

18 

Computed  semi-minor  axis  length: 

SMI=val. 

CONT 

19 

Computed  containment  ellipse  area: 
See   Note  4. 

A  =  val. 

20 

For  a  containment  ellipse  of  a  given 
size,  first  press  Fl,  0. 

Fl  0 

21 

Next,  key  in  4  and  press  EXE. 

OPTION? 

4  EXE 

22 

Key  in  the  desired  containment  ellipse 
size  (see  Note  3) . 

SIZE? 

EXE 

23 

Computed  value  of  the  containment 
probability: 

PRB=val. 

CONT 

24 

Computed  semi-major  axis  length: 

SMJ=val. 

CONT 

25 

Computed  major  axis  direction: 

DRC=val. 

CONT 

26 

Computed  semi-minor  axis  length: 

SMI=val. 

CONT 

27 

Computed  containment  ellipse  area: 
See   Note  4. 

A  =  val. 

28 

To  enter  additional  data  from  new  or 
old  stations,  first  press  Fl,  0. 

Fl  0 

29 

Next,  key  in  5,  press  EXE  and  then 
repeat  Steps  3,  4,  5  &  6. 

OPTION? 

5  EXE 

11 


Notes : 

1.  Enter  the  program  in  P0  for  Fl  0  activation.    Before  running 
the  program,  first  press  F2  then  DEFM  then  1  and  then  EXE. 
This  is  required  in  order  to  use  the  array  variables  in  the 
program. 

2.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

3.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 

probability  of  containment,  the  minimum  area  containment 

region  is  an  ellipse  centered  on  the  position  estimate.   The 

semi-major  axis  =  ka,,^.  and  semi-minor  axis  =  ka,..,.  where  k 
J  MJ  MI 

is  the  size  of  the  ellipse  and  a.,_  and  o.tT    are  the  standard 

c  MJ       MI 

deviations  (uncertainty  measure)  of  the  position  estimate 
in  the  major  axis  and  minor  axis  directions. 

4.  The  area  units  are  the  range  units  squared. 

5.  For  a  definition  of  the  program  options  and  their  function, 
press  Fl,  0.   Then  after  OPTION?  is  displayed,  enter  0  and 
press  EXE.   The  display  will  show:   TPE  =  1  EST  =  2  SIZ  =  3. 
Next,  press  CONT.   The  display  will  show:   PRB  =  4  CON  =  5. 
To  repeat  the  displays,  press  CONT. 


12 


POSITION  ESTIMATION  PROGRAM 
TI-59  USER  INSTRUCTIONS 


Step 

Instruction 

Enter 

Press 

Display 

1 

Enter  the  program  (see  Note  1) . 

2 

To  run  the  program,  press  A. 

A 

9 

3 

Enter,  in  decimal  degrees,  the  observed 
bearing  of  the  target  from  a  station  or 
the  reciprocal  of  the  observed  bearing 
of  the  station  from  the  target  (see 
Note  2) . 

0 

B 

e 

4 

Enter  the  station  bearing  in  decimal 
degrees  from  the  reference  point.   Use 
zero  if  the  station  is  the  reference 
point. 

a 

R/S 

a 

5 

Enter  the  station  range  from  the  refer- 
ence point.   Use  zero  if  the  station  is 
the  reference  point.   Use  any  units. 

P 

R/S 

P 

6 

Enter  the  bearing  error  (standard  devi- 
ation) in  decimal  degrees.   After  press- 
ing R/S ,  the  display  indicates  the  order 
number  of  the  data  entry. 

e 

R/S 

n 

7 

Repeat  Steps  3,  4,  5  &  6  for  at  least 
one  more  station. 

8 

Display  a  bearing  estimate  $   with  respect 
to  the  reference  station. 

C 

/s 

9 

Display  a  range  estimate  r  with  respect 
to  the  reference  station.   For  an  ellip- 
tical containment  region,  go  to  Step  10 
or  Step  15.   To  enter  additional  data 
from  new  or  old  stations,  repeat  Steps  3, 
4,  5  &  6. 

R/S 

r 

.0 

For  a  containment  ellipse  of  a  given  con- 
tainment probability,  enter  the  contain- 
ment probability.   Next,  press  E  and 
display  the  ellipse  size  (see  Note  3). 

P 

E 

k 

.1 

Display  the  semi-major  axis  length. 

R/S 

MJ 

.2 

Display  the  major  axis  direction. 

R/S 

y 

.3 

Display  the  semi-minor  axis  length. 

R/S 

kaM, 
MI 

13 


Step 


Instruction 


Enter 


Press 


Displ 


Display  the  containment  ellipse  area. 
(See  Note  4.) 

For  a  containment  ellipse  of  a  given 
size,  enter  the  containment  ellipse 
size  (See  Note  3) .   Next,  press  D  and 
display  the  containment  probability. 

16     Display  the  semi-major  axis  length. 

Display  the  major  axis  direction. 

Display  the  semi-minor  axis  length. 

Display  the  containment  ellipse  area. 
(See  Note  4.) 

To  enter  additional  data  from  new  or 
old  stations,  repeat  Steps  3,  4,  5 
&  6. 


R/S 


area 


R/S 
R/S 
R/S 
R/S 


k0MT 

MJ 
Y 

kaMIj 
area 
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Notes: 

1.  The  program  requires  the  normal  partition.   If  the  calculator 
has  been  in  use,  this  can  be  assured  by  turning  the  calculator 
off  and  then  on  before  loading  the  program. 

2.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

3.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 
probability  of  containment,  the  minimum  area  containment  region 
is  an  ellipse  centered  on  the  position  estimate.   The  semi- 
major  axis  =  ko,.T  and  semi-minor  axis  =  ka.,_  where  k  is  the 

.  J  MJ  MI 

size  of  the  ellipse  and  o,,^.   and  a.,_  are  the  standard  deviations 

^  MJ       MI 

(uncertainty  measure)  of  the  position  estimate  in  the  major 
axis  and  minor  axis  directions. 

4.  The  area  units  are  the  range  units  squared. 

5.  Negative  bearing  estimates  and  negative  major  axis  directions 
can  result.   To  convert  a  negative  bearing  estimate  to  the 
value  that  would  be  output  by  the  other  programs,  add  360°  to 
the  estimate.   For  example,  -5°  becomes  355°.   To  convert  a 
negative  direction,  add  180°.   For  example,  -5°  becomes  175°. 
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III.   Two  Examples 

In  Scenario  1,  the  scenario  for  the  first  example,  bearings 
are  taken  on  a  target  from  three  separate  stations.   Figure  1  on 
Page  18  shows  Scenario  1  and  Table  1  below  gives  the  station  data. 
The  stations  are  numbered  according  to  the  order  of  station  data 
input  to  the  program. 


TABLE  1 

OBS  BRG 

STA  BRG 

STA  RNG 

BRG 

ERR 

Station 

1 

038° 

334° 

13500 

4° 

Station 

2 

324° 

050° 

11350 

3° 

Station 

3 

003° 

000° 

0 

4° 

Note  from  Table  1  that  the  reference  point  is  at  Station  3. 

Program  outputs  for  Scenario  1  are  given  in  List  1  on  Page 
19.   List  1  and  List  2  which  gives  program  outputs  for  Scenario 
2,  are  copies  of  printer  tapes  that  were  generated  with  a  Casio 
FP-10  printer  and  a  Casio  FX-702P  using  the  Casio  FX-702P  program. 
With  allowance  for  differences  in  user  instructions,  display  for- 
mats and  round-off  errors,  the  tapes  indicate  the  output  that 
should  be  obtained  using  any  of  the  other  calculator  programs  ex- 
cept for  the  TI-59  program  where  equivalent  negative  angles  occur. 

In  Scenario  2,  the  scenario  for  the  second  example,  bearings 
are  taken  from  a  target  on  three  separate  stations.   Figure  2  on 
Page  20  shows  Scenario  2  and  Table  2  below  gives  the  station  data. 
As  in  Scenario  1,  the  stations  are  numbered  according  to  the  order 
of  station  data  input  to  the  program.   In  this  scenario  the  reci- 
procal of  the  observed  bearings  are  used  in  order  to  provide  an 
equvalent  scenario  that  is  appropriate  for  the  program. 
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TABLE  2 

OBS  BRG     RCP  BRG     STA  BRG  STA  RNG  BRG  ERR 

Station  1        211°        031°        000°          0  3° 

Station  2        172°        352°        115°  11100  3° 

Station  3        146°        326°        082°  13800  3° 

Note  from  Table  2  that  the  reference  point  is  at  Station  1. 

Program  outputs  for  Scenario  2  are  given  in  List  2  on  Page 
21.  The  two  bearing  estimates  in  the  list  are  the  reciprocal  of 
the  target's  estimate  of  the  bearing  of  Station  1. 


17 


Position  Estimate 


Station  2 


Station  3 


Figure  1.   A  position  estimate  and  a  .33  confidence  region. 

Table  1  gives  the  bearings  of  the  indicated  bearing 
lines  and  List  1  gives  the  bearing  and  the  range  of 
the  position  estimate  from  Station  3. 
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LIST  1 

OPTION? 

OPTION? 

TPE=1  EST=2  SI2=3 

OBS  BRS? 

7 

PRB=4  C0H=5 

STfl  BRS? 

OPTION? 

1 

OBS  8R6? 

TO 

STfl  BR6? 
334 

STfl  RN6? 
13599 
BR6  ERR? 

STfl  RNS? 

6 

BRS  ERR? 

4 

OBS  8R8? 

OPTION? 

2 

BR6=  9.15 

OBS  3R6? 

324 

STfl  BR6? 

59 

STH  RH6? 

11359 

BR6  ERR? 

RH6=  19553.78 

END 

OPTION? 

7 

SIZE? 

OBS  8R6? 

PRB=  9.86 

OPTION? 

SRJ=  1712.95 

BR6=  359.51 

DIR=  12.48 

RN8=  19494.39 

SHI=  1129.99 

END 

fi=  6889418.49 

OPTION? 

END 

SIZE? 
PRB=  9.86 

OPTION? 

4 

PRB? 

SHJ=  1737.32 

SIZE=  2.15 

DIR»  17.69 

SHJ=  1837.97 

SHI=  1232.96 

DIR=  12.48 

H=  6729444.91 

Sfil=  1212.36 

END 

fi=  7998348.38 

OPTION? 

4 

PRB? 

.9 

SIZE=  2.15 

END 

SHJ=  1364.12 

DIR=  17.69 

SHI=  1322.94 

H=  7747559.77 

END 
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Position  Estimate 


Target 


Station  3 


Station  2 


Figure  2.   A  position  estimate  and  .71  confidence  region. 

Table  2  gives  the  bearings  of  the  indicated  bearing 
lines  and  List  2  gives  the  reciprocal  bearing  and 
the  range  of  the  position  estimate  from  Station  1. 
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List  2 


OPTION? 

ft 

TPE=1  EST=2  SIZ=3 
PRBM  C0N=5 

OPTION? 

t 

OBS  BR6? 

31 

STfi  BRG? 

H 

STfi  RHG? 

8 

BR6  ERR? 

7 

OBS  BRG? 

352 

STfi  BRG? 

115 

STfi  RHG? 

11189 

BRG  ERR? 

3 

OBS  BRG? 

OPTION? 

2 

BRG=  31.88 

RN8=  14792.53 

END 

OPTION? 
3 

SIZE? 

1 

PRB=  8.86 

S«J=  3686.56 

DIR=  15.39 

SHI=  1253.72 

fl=  14285978.91 

END 

OPTION? 
4 

PRB? 
,9 

SIZE=  2.15 
SHJ=  3369,7? 
DIR=  15.39 
SHI=  1345.22 
fl=  16354192.26 


OPTION? 

r 

OBS  BRG? 

326 

STfi  BRG? 

82 

STfi  RHG? 

13888 

BRG  ERR? 

3 

OBS  BRG? 

OPTION? 
•i 

BRS=  31.65 
RN8=  13765.38 

END 
OPTION? 

7 
■i 

SIZE? 

2 

PRB=  8.86 

SHJ=  2323.19 

DIR=  172.86 

SHI=  967.98 

fi=  6151933.32 

END 

OPTION? 

4 

PRB? 

3!ZE=  2.15 

SHJ=  2178.85 

DIR=  172.86 

Sftl=  1338.54 

R=  7882732.54 

END 


IV.       Program  Listings 
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HP-41CV    PROGRAM 


81*LBL  "TPE" 
92  DEG 

03  CLRG 

04  .801 
85  STO  89 
06  26 

87  STO  18 
83*LBL  "COH" 
09  FIX  0 

18  -DBS  BRG?" 

11  PROMPT 

12  STO  11 

13  "STfl  BRG?" 

14  PROMPT 

15  STO  12 

16  -STfi  RNG?" 

17  PROMPT 
13  STO  13 

19  "BRG  ERR?" 

28  PROMPT 

21  STO  14 

22  1 

23  ST+  88 

24  2 

25  RCL  83 

26  X>Y? 

27  GTO  82 
23  1 

29  ST+  18 

38  RCL  11 

31  STO  IHD  18 

32  1 

33  ST+  18 

34  RCL  12 

35  STO  IHD  10 

36  1 

37  ST+  18 
33  RCL  13 

39  STO  IND  10 
48  1 

41  ST+  18 

42  RCL  14 

43  STO  IND  18 

44  1 

45  RCL  88 

46  X=Y? 

47  GTO  "COH" 

48  RCL  29 

49  RCL  28 
58  RCL  27 


51  - 


52  SIN 

53  * 

54  STO 

15 

55  RCL 

33 

56  RCL 

57  RCL 

31 

58  - 

59  SIN 

68  * 

61  STO 

16 

62  RCL 

71 

63  RCL 

97 

64  - 

65  SIN 

66  STO 

17 

67  X=8? 

63  GTO 

88 

69  RCL 

15 

78  RCL 

31 

71  SIN 

72  * 

73  RCL 

16 

74  RCL 

C  i 

75  SIN 

76  * 

77  - 

78  RCL 

17 

79  / 

38  STO  81 

81  RCL 

15 

32  RCL 

31 

33  COS 

84  * 

85  RCL 

16 

36  RCL 

»37 

L  i 

37  COS 

33  * 

89  - 

98  RCL 

17 

91  / 

92  STO 

ez 

93  26 

94  STO 

18 

95*LBL  85 

96  1 

97  ST+ 

18 

93  RCL 

IND  18 

99  STO 

11 

188  1 

101 

ST+ 

18 

102 

RCL 

IND  18 

183 

STO 

12 

184 

i 

i. 

105 

ST+ 

18 

106 

RCL 

IND  18 

187 

STO 

13 

183 

1 

109 

ST+ 

18 

110  RCL 

IND  18 

111 

STO 

14 

112 

XEQ 

87 

113 

ISG 

89 

114 

GTO 

85 

115 

GTO 

•CON" 

lifo 

LBL 

82 

117 

XEQ 

67 

113 

GTO 

"CON" 

H9*LBL 

9?" 

128  RCL 

81 

121 

RCL 

13 

122 

RCL 

12 

123 

SIN 

124 

* 

125 

- 

126 

RCL 

92 

127 

RCL 

|  7 

i  --1 

123 

RCL 

12 

129 

COS 

138 

* 

131 

- 

132 

R-P 

133 

STO 

16 

134 

y<  >y 

135 

STO 

15 

136 

RCL 

11 

137 

RCL 

(  c 

i  J 

138 

- 

139 

138 

148 

X<>' 

' 

141 

,V\  - 

n 

142 

GTO 

u 

143 

368 

144 

- 

145«LBL 

ii 

146 

RCL 

14 

147 

148 

STO 

1 1 

149  RCL 

16 

158  RCL 

14 

23 


151  * 

152  PI 

153  * 

154  188 

155  / 

156  STO  18 

157  RCL  15 

158  COS 

159  RCL  13 

168  / 

161  STO  19 

162  RCL  15 

163  SIN 

164  RCL  18 

165  / 

166  STO  28 

167  RCL  19 

163  Xt2 

169  ST+  83 

178  RCL  19 

171  RCL  28 

172  * 

173  ST+  84 

174  RCL  28 

175  Xt2 

176  ST+  85 

177  RCL  17 

178  RCL  19 

179  * 

188  ST+  86 

181  RCL  17 

132  RCL  28 

133  * 

184  ST+  87 

135  RTN 

136*LBL  "ES 

187  FIX  2 

188  RCL  84 

139  Xt2 

190  RCL  03 

191  RCL  85 

192  * 

193  - 

194  STO  28 

195  X=8? 

196  GTO  88 

197  RCL  81 

198  RCL  84 

199  RCL  87 

200  * 

201  RCL  05  251  COS 

202  RCL  06  252  Xt2 

283  *  253  RCL  85 

284  -  254  * 

285  RCL  28  255  RCL  23 

286  /  256  1 

287  +  .257  P-R 
283  STO  26  253  * 
289  RCL  82  259  2 

218  RCL  83  268  * 

211  RCL  87  261  RCL  04 

212  *  262  * 

213  RCL  04  263  STO  26 

214  RCL  06  264  - 

215  *  265  RCL  23 

216  -  266  SIH 

217  RCL  28  267  Xt2 
213  /  268  RCL  83 

219  +  269  * 
228  RCL  26  278  + 

221  XOY  271  RCL  28 

222  R-P  272  CHS 

223  STO  21  273  / 

224  XOY  274  STO  24 

225  X>8?  275  RCL  23 

226  GTO  86  276  SIN 

227  368  277  Xt2 
223  +  273  RCL  85 
229»LBL  86  279  * 

238  STO  22  238  RCL  26 

231  RCL  84  281  + 

232  RCL  83  282  RCL  23 

233  RCL  85  233  COS 

234  X*Y?  234  Xt2 

235  GTO  83  285  RCL  83 

236  RCL  84  286  * 

237  SIGN  237  + 

238  45  283  RCL  28 

239  *  289  CHS 
248  GTO  84  298  / 
24ULBL  83  291  STO  25 

242  -  292  RCL  24 

243  /  293  X<=Y? 

244  2  294  GTO  81 

245  *  295  STO  25 

246  ATAN  296  XOY 

247  2  297  STO  24 

248  /  298  98 
249*LBL  84  299  ST+  23 
258  STO  23  388*LBL  81 
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381  "BRG=- 

382  ARCL  22 

383  AVI EH 

384  STOP 

385  "RNG=" 

306  flRCL  21 

307  AVIEH 

308  STOP 

309  'END- 

310  GTO  08 
311*LBL  -SI2- 

312  "SIZE" 

313  PROMPT 

314  STO  11 

315  Xt2 

316  2 

317  / 
313  CHS 
319  EtX 
328  1 

321  XOY 

322  - 

323  STO  88 

324  -PRB=" 

325  flRCL  80 

326  AVI EH 

327  STOP 

328  GTO  89 
329*LBL  *PRB" 
338  "PRB?" 

331  PROMPT 

332  CHS 

777  t 

-jjj  1 

334  + 

335  LH 

336  2 

337  * 

338  CHS 

339  SQRT 
348  STO  11 

341  -SIZE=- 

342  ARCL  11 

343  AVI EH 

344  STOP 
345»LBL  89 

346  RCL  23 

347  STO  12 

348  RCL  11 

349  RCL  25 

350  SQRT 


351 

* 

352 

RCL  11 

353 

RCL  24 

354  SQRT 

355 

* 

356 

XOY 

357 

-SHJ=" 

358 

ARCL  X 

359 

AVIEH 

368 

STOP 

361 

RCL  23 

362 

X>8? 

363 

GTO  18 

364 

138 

365 

+ 

366*LBL  18 

367 

■DRC=" 

368 

ARCL  X 

369 

AVIEH 

378 

STOP 

371 

-SHI=" 

372  RDN 

373 

flRCL  Y 

374  AVI EH 

375 

STOP 

376 

PI 

777 

* 

373 

* 

379 

■A=" 

388  ARCL  X 

381' 

HBL  88 

382  AVIEH 

333 

STOP 

384 

-END" 

335  GTO  88 

386*LBL  88 

387 

AVIEH 

388 

STOP 

389 

-HO  SOL 

398  GTO  83 

391 

END 
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SHARP  PC-1500  PROGRAM 


18:  "A":CLEAR  :DIP1 
A (7): DEGREE 

15: "C": INPUT  "OBS 
BRG?  ";P: 
INPUT  "STA  BRG 
?  ",'Q: INPUT  "S 
TA  RNG?  "jR: 
INPUT  "BRG  ERR 
?  "JO: PAUSE  " 

20: IF  I=2GOTO  50 

25: I=I+1:AC1-1)-P 
:A(I+1)=Q: A(I  + 
3)=R:A(I+5)=0: 
IF  I^IGOTO  15 

30:X=A(4)*SIN  (AC 
2)-A(0)): Y=A(5 
)*SIN  (A(3)-A( 
1)):Z=S1N  (A(l 
)-A(0)): IF  Z=0 
GOTO  150 

35:U=(X*SIN  A(l)- 
Y*SIN  A(0))/Z: 
U=(X*COS  A(l)- 
Y*COS  fi(0))/Z 

40: FOR  n=0TO  1 

45:P=ACn):Q=A(n+2 
):R=A(n+4):0=A 
(M+6):G0SUB  12 
0:NEXT  n:GOTO 
15 

50:GOSUB  170: GOTO 
15 

60: "Z": PAUSE  "  " : 
F=(B*B-A*C): IF 
F=0GOTO  150 

65:X=U+(B*E-C*DV 
F: Y=U-KA*E-B*D 
)/F:GOSUB  200 

70: T=SGN  B*45: IF 
A=CGOTO  80 

75:T=.5*ATN  (2*B/ 

ca-o; 

80:g=(c*cos  t*cos 

T-2*BXC0S  T* 
SIN  T+A*S1N  T* 
SIN  T)/-F:G=TG 

85:H=(C*SIN  T*SIN 
T+2*B*C0S  T* 
SIN  T+A*COS  T* 
COS  T)/-F:H^TH 
: IF  H>^GGOTO  9 
5 

90: Z=H:H=G:G=Z: T  = 
T+90 


95:PRINT  "BRG=";J 

.'PRINT  "RNG="; 

K:G0T0  145 
100: "H":PRINT  "  TP 

E=A  EST^Z  SIZ= 

S  PRB=X" 
105: PRINT  "CON^C": 

GOTO  100 
120: "S": INPUT  "SIZ 

E?  ";S:0^1-EXP 

<-S*S/2) 
125:PRINT  USING  "tt 

#.«**";  "PRB=";0 

:US1NG  :GOTO  1 

35 
130: "X": INPUT  "PRB 

?  ";0:S=T(-2* 

LN  (1-0)): 

PRINT  "SIZE^"; 

S 
135:X=S*G:Y=S*H:N= 

T:PR1NT  "SriJ=" 

JY: IF  N<0LET  N 

=N+180 
140:PRINT  "D1R=";N 

:PR1NT  "sni="; 

X:PRINT  "A=";n 

*X*Y 
145:PRINT  "END": 

GOTO  145 
150:PRINT  "NO  SOL" 

:GOT0  150 
170:X=U-R*S1N  Q:Y^r 

U-R*COS  Q: 

GOSUB  200 
175:U=(P-J):L^K*0* 

JT/180:G-COS    J/ 

l_:H=SIN    J/L:  IF 

W>=180LET    U=U- 

360: GOTO  185 
180: IF  U<=-180LET 

U=UI+360 
185:U=W/0: A=G*G+A: 

B=G*H+B:C=H*H+ 

C:D=UIXG+D:  E=W* 

H+E: RETURN 
200:K=T(X*X+Y*Y) : 

IF    K^0LET    J=0: 

RETURN 
205: J^ACS    (Y/K): IF 

ASN     (X/KX0LET 

J=360-J 
210: RETURN 
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SHARP    PC-1211    PROGRAM 


10: "A": CLEAR   : 
DEGREE 

15: "C": INPUT  "0 
BS  BRG?  "5P: 
INPUT  "STA  B 
Fij''  "SQ: 
INPUT  "STA  R 
NG?  "5R: 
INPUT  "BRG  E 
FR?  "50: IF  1 
=2G0TCi   130 

75: I=I+liA<I+26 
>=P: A<I+28)= 
Q:Aa+30>=R: 
A<I+32>=0:IF 
1=1 GOTO    15 

85:K=A<31>+SIN 
CA<2$0-A<:27> 
):Y=A<32>* 
SIN    CA-30)-A 
C28>>:Z=SIN 
CA(28)-AC27) 

'30:  IF  2=OG0TCi  3 

25 
95:  U=*<X*SIN   A<2 

8>-Y*SIN   A<2 

COS   A<28>-Y* 

COS  A<27>)/2 
105: FOR  I=1T0  2: 

F*=A<I+26>:Q= 

A<I+28>:R=A< 

I+30>:0=A<I+ 

32): GOSUB  40 

0 
110: NEXT    I: GOTO 

15 
130:G0SUB  400 
135: GOTO    15 
140:  "Z":F=<B*B-A 

*C>: IF  F=0 

GOTO   325 
150:X=U+CB+E--C+B 

>,'F:V=V-KA*E 

-B*D>-'F: 

GOSUB  500 
160:T=SGN  B*45: 

IF  A----CTHEN   1 

70 
165:T=.5*ATN    (2* 

B-'<A-C>> 

170:  6=<C*C05    T* 
COS   T-2*B* 
COS   T+SIN  T+ 
A*SIN  T+SIN 
T>.'-F:G=TG 

175:H*<C*SIN  T* 
SIN  T+2-+B* 
COS  T*SIM  T+ 
A*COS  T*CGS 
T;.-'-F:  H=TH: 
IF  H>=GGOTO 
135 


180:Z=H:H--G-*G=Zs 
T-T+90 

135: PRINT  "ERG=" 
5 J : PRINT  "RN 
G:="sK:G0T0  3 
20 

200: "H": PRINT  "T 
PE=A  EST=Z  S 
IZ=S   PRB--X" 

205: PRINT  "CON*C 
": GO TO  200 

225:  "•:-:.":  INPUT  "S 
IZE?  "5S:0=1 
-EXP    <-S*S--2 

230: PRINT   USING 
"*«. ttlt" ; "PRB 
="5  0: USING   : 
GOTO  300 

235: "X": INPUT  "P 
RB?  "J0:5=r< 
-2+LN  <l-0>) 
: PR I NT  "SIZE 
-"jS 

30Q:X=S*G:Y=S*H: 
N=T: PRINT   "S 
MT=" ;  Y:  IF  N 
OOLET  N--N+1 
80 

315: PRINT  "DRi>" 
SN: PRINT  "SN 
I=,,JX:Z=Jt*X* 
V:  PR  INT  "A=H 

■  T 

320: PRINT  "END": 
GOTO  320 

325: PRINT  "NO  SO 
L":GOTO  325 

40Q:X=U-R*SIN  0.: 
V= V-  R+CO  S  Q : 
GOSUB  500: Ws 
(P-J>:L-K*0* 
jr,-i30:G:-C0S 
J.-L:H=SIN  J/ 
L 

415: IF  U>=130L£7 
W=W-36Q: GOTO 
440 

425:  IF  NO- 180 
LET   W*W+360 

440:  W=W--'0:A=G*G+ 
h:B-G+H-+B:C= 
H*H+C:B=W*G+ 
D:E=H*H+E: 
RETURN 

500 :  K=T  ■■  X*X+  Y*Y ) 
: IF  K=OLET  J 
=0: RETURN 

510:  M- ASH  CX/'K): 
J=ACS  CY.-'K): 
IF  hKOLET  J= 
360- J 

515:  RETURN 
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CASIO    FX-702P    PROGRAM 


5  IHP  "OPTION", Ji  i25  H=(C*SIN  T*SIH 

SET  F2:H0DE  4  T+2*B*C0S  T*SIH 

18  IF  J=8  THEN  175  T+fl*C0S  T*C0S 

15  IF  J=2  THEN  188  T)/-F 

28  IF  J=3  THEN  145  138  8=SQR  G;H=SQR  H 

25  IF  J=4  THEN  158  : IF  S>H?Z=H:H=6 

38  IF  J=5  THEN  45  :G=Z:T=T+90 

48  VflC  :SET  F2  135  IF  Y<8;Y=Y+3b8 

45  IHP  "DBS  BR6",P  148  PRT  "BRG=";Y;PR 

J  INP  ??STfi  BRG",  T  3RHG=5';*:G0T0 

Q: IHP  "STfl  RN6"  165 

,R  145  IHP  nSI2E5',S:0= 

58  INP  "BR6  ERRw,0  i-EXP  f-S*S/2): 

:IF  1=2  THEN  95  PRT  "PPR=";!i:fif! 

55  I=I+i:fl(I-l)=P;  TO  155 

fi(I+l)=Q:fl(I+3)  158  INP  WPRBS* 0: S=S 

=R:fl(I+5)=0  §R  (-2*LH  (1-0) 

68  IF  1=1  THEN  45  ):PRT  "SIZE=B;S 

65  X=fi(4)*SIN  (fi(2  155  ;<=S*G:Y=S*H:FRT 

)-fi(8)):Y=H(5)*  "SHJ=7f;  Y:N="t  I 

SIN  (fi(3H(i>)  F  H<8;H=Nti86 

78  Z=SIN  (H(1H(8  168  PRT  nDIR=B;NiPR 

)>:IF  2=8  THEN  T  "SHI^WPRT 

170  ffR=w;ji*X*Y 

75  U=(X*SIH  H(l)-Y  165  PRT  "END*:GOTO 

♦SIN  fl(8»/2:V=  165 

(X*COS  R(1M*C  178  PRT  "NO  SOL3:  GO 

OS  fi(8))/Z  TO  176 

88  FOR  H=0  TO  l:P=  i75  PRT  "TPE=1  EST= 

fi<M):Q=fl(M+2)  2  SIZ=3" 

35  R=fi(«+4):0=fl(H+  188  PRT  "PRB=4  CON= 

fa-sGSB  288: NEXT  5'J:S0T0  175 

H:aOT0  45  288  X=U-R*SIN  S:Y=V 

98  R=H(M+4):0=fi(H+  -R*C03  Q:RPC  Y, 

6):GSB  288: NEXT  I 

mm  45  285  H=(P-Y):L=X*0*s 

95  SSB  288: GOTO  45  /186:G=C0S  Vl\ 

188  F=(B*6-H*C):IF  H=SIN  VI 

F=8  THEN  178  210  IF  Hil88iH=H-36 

105  X=U+<B*E-C*D)/F  8 : GOTO  220 

:Y=V+(fl*E-B*D)/  215  IF  !te-188;Ni+3 

F:RPC  Y,X  68 

118  T=SGN  B*45:IF  H  228  H=H/0:R=6*G+fiiB 

=C  THEN  128  =6*H+B:C=H*H+C: 

115  T=.5*fiTN  (2*B/(  D=w*G+D:E=i»i*H+E 

fl-0)  :RET 
120  6=(C+C0S  T*C0S 
T-2*B*C0S  T*9IN 
T+fl*SIN  T*SIN 
T)/-F 
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Appendix  1.   The  Estimation  Procedure 

In  this  development  the  assumptions  stated  in  Section  I  are 
required  conditions.   A  rectangular  coordinate  system  is  used 
with  the  positive  y-axis  directed  north,  the  positive  x-axis 
directed  east  and  the  origin  at  the  reference  point  and  all  angles 
are  in  radians.   Figure  3  shows  a  station  located  with  respect  to 
the  coordinate  system.   The  bearing  line  of  length  r  goes  to  the 
object's  unknown  position,  the  bearing  line  of  length  r  goes  to 
an  initial  estimate  of  the  object's  position  and  the  third  bear- 
ing line  corresponds  to  an  observed  bearing. 

To  determine  the  coordinates  for  a  final  estimate,  consider 
the  arc  segments  u  =  r(0-(f))  where  9-$  is  the  bearing  error  and 
v  =  r(cj>-3)  and  w  =  r(9-3)  that  are  defined  by  the  three  bearing 
lines  and  the  circle  of  radius  r  that  is  centered  on  the  station 
and  goes  through  the  initial  estimate.   The  geometry  is  shown  in 
Figure  3.   Note  that  u  can  be  defined  by  u  =  w  -  v .   In  this 
expression,  w  =  r(9-3)  is  known,  and  v  can  be  determined  in  terms 
of  x  and  y  the  unknown  coordinates  of  the  target.   With  the  refer- 
ence point  at  the  initial  estimte,  to  first  order, 
v  =  x  cos  3  -  y  sin  3  and  u  =  w  -  x  cos  3  +  y  sin  3.   Since  this 
approximation  applies  to  all  stations,  its  use  suggests  that, 
for  each  station,  r  -  r  which  is  equivalent  to  having  the  initial 
estimate  relatively  close  to  the  target's  position. 

Since,  for  each  station  i,  an  observed  bearing  9^  is  the 

value  of  a  normal  random  variable  0.  with  mean  <j)^f  the  coordinate 

u.  =   r.  (9. -<}).)  is  the  value  of  a  normal  random  variable  U.  with 
l      11   l  1 

mean  zero.   In  addition,  since  the  0^^  are  independent,  the  IT  are 
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Figure  3.   The  coordinate  geometry.   The  coordinates  of  the  initial 

estimate  are  (x*,y*).   In  the  development,  the  reference 

point  is  at  the  initial  estimate. 


34 


also  independent.   And,  consequently,  the  joint  distribution  of 
the  U.  is  determined. 

To  estimate  x  and  y,  maximum  likelihood  estimates  are  used 
here.   The  likelihood  for  a  sample  Q^,    62,...,  9   from  n  stations 


is 


n     1  1   n         2   ? 

Me,  ,e0,.  ..,8)  =        n  — - —   exp  -  ^  z  (e.-<j>. )  /e 

12      n       1  /2?  e.  2   1    x   x    x 


and  the  likelihood  for  a  cooresponding  sample  u, ,  u~,...,  u   is 

n  1   n   2   2 

L  (u,  ,u„ ,  .  .  .  ,u  )  =    IT  exp  -  ■*•  Z  u./a. 

1  /2tt  a.  1 

i 

with  a.  =  r.e.  where  e.  is  the  standard  deviation  of  0 . . 
1111  l 

By  definition,  the  maximum  likelihood  estimates  for  x  and  y 

are  the  estimates  which  make  L(u, ,  u^,...,u  )  a  maximum.   In 

1    2       n 

this  case,  making  L (u, ,u? , . . . ,u  )  a  maximum  is  equivalent  to 

n        2        2 

making      Z(u.    /a.    )    a  minimum.      So,    to    find   the   maximum   likelihood 

1    ^        \ 
estimates   x   and   y,    solve    the    following   two   equations    for   x   and   y: 

3(*n  L)|      =   0        and        ii£LLL|      =   o 

9x         '  3y         ' 

x   =   x  y    =   y 


With   u.    =  w.    -   x   cos    3.    +   y    sin    p.,    the   equations   can   be   written 


as    follows: 


and 


n  2 

Z    (w.    -   x  cos    3.    +   y   sin   3.)  (cos    3^/ck    =   0 


n  2 

Z    (w.    -   x   cos    3-    +   y   sin    3-) (sin    3-)/a.    =    0 
l  i  i  ii 
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In  terms  of  the  following  quantities: 


2       2  2 

A  =  Z(cos   $.)/a.  ,  B  =  I (sin  3-  cos  3 • ) /a   , 


C  =  I  (sin2  (3.) /a2  ,  D  =  Z(w.  cos  3.) /a?  , 

i   i  i      i  '    i  ' 


E  =  I (w.  sin  3 ■ )/a2  , 
i      li 


the  equations  become: 

Ax  -  By  =  D 


s\  /\ 


Bx  -  Cy  =  E 
The  solutions  are: 


^ 


2 


(1)  x  =  (BE  -  CD)/(B   -  AC) 
and 

(2)  y  =  (AE  -  BD)/(B2  -  AC) 

A  confidence  region  can  be  constructed  about  an  estimated 
position.   In  order  to  indicate  how  this  can  be  done,  a  proba- 
bility region  about  the  true  position  will  be  considered  first 

Both  x  and  y  are  values  of  random  variables.   If  a  new  set 
of  bearings  0,,  92/  •••/  6   is  observed  (for  the  same  initial 
estimate  and  a  fixed  target) ,  in  general,  a  new  pair  of  values 
x  and  y  will  be  obtained. 
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If  X  and  Y  represent  these  random  variables, 


1 

X  = 


n 


5 E  (W./a.)(B  sin  0.  -  C  cos  3-) 

1 -AC)   \         ^      x  x  ! 


(B  -AC) 


1 

Y  = 


n 


~ Z  (W./o  ) (A  sin  3-  -  B  cos  3-) 

-AC)   1    1   1  x  ± 


(B  -AC) 


where  W.  =  r.(0.  -  3.). 
l     11     i 

Since  X  and  Y  are  a  linear  combination  of  the  n  normal  random 
variables  V^,  W2 ,  ...,  W  ,  or  equivalently  of  the  n  normal  random 
variables  0^,    02,  ...,  o  ,  they  have  a  joint  normal  distribution. 
Since  E(Wi)  =  ^((Jk  -  3±)  ,   if  6j_  =  ^  for  i  =  1,  2,  ...,  n,  that 
is,  if  the  initial  estimate  of  the  target's  position  is  at  the 
target's  position,  E(W1)  =  0  for  i  =  1,  2,  ...,  n.   In  this  case 
E(X)  =0  and  E (Y)  =0  and  the  joint  normal  distribution  is  centered 
on  the  object's  position.   To  the  degree  of  the  approximations 
that  have  been  made,  this  is  also  true  if  the  initial  estimate  is 
not  at  the  target's  position. 

A  region  of  minimum  area  for  a  given  probability  of  contain- 
ment of  an  estimated  position  can  be  determined.   The  region  is 
bounded  by  an  ellipse  which  is  centered  on  the  object's  position 
and  whose  axes  lie  along  the  axes  of  an  x'y' -coordinate  system 
that  is  obtained  by  rotating  the  xy-coordinate  system  that  is 
centered  on  the  object's  position  through  an  angle  y.      In  this 

system,  a^,  A,  is  0,  that  is  x'  and  y'  are  independent  normal  ran- 
x  y 

dom  variables.   The  two  coordinate  systems  are  illustrated  in 
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Figure  4.   The  coordinates  of  a  point  in  the  two  systems  are 
related  by 


x'  =  x  cos  y  -  y  sin  y 


y1  =  x  sin  Y  +  y  cos  Y 


These  relations  imply 


2      2     2  2     2 

(3)    o^   =  a„  cos   y  -  2a^  cos  y  sin  y   +   a^  sin  y    , 

xy  y 


2      2     2  2     2 

(4)    a^   =  a^  sin  y   +   2a^  cos  y   sin  y   +   o^   cos  y 

y1     x  xy  y 


and 


2         2     2  2         2 

(5)    a^    A   =  (aA  -  a  J  sin  y   cos  y  +  a^  (cos  y    -    sin   y) 

x'y1      x    y  xy 


where  y,  the  angle  of  rotation  of  the  coordinate  axes,  is  posi- 
tive in  the  clockwise  direction.   And  o",*.    =    0  implies 

x'y1        r 

2oaa 
tan  2y  = 


2     2 

o      -    o 


y   x 

With  the  initial  estimate  of  the  target's  position  at  the 

2 
target's  position  E(W.)  =0  and  therefore,  Var (W .)  =  o.       for 

i  =  1,  2,  ...,  n.   In  this  case 

2       1      n     2  2 

o^    =  — ^ =■   Z  (l/o.)   (B  sin  3-  -  C  cos  B. )  , 

x    (B  -AC^   1      1  1 
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North 


Reference 
Point 


East 


Figure  4.   An  elliptical  confidence  region  and  the  primed 
coordinate  system  in  which  the  covariance  ^,\ 
is  zero.   The  center  of  the  ellipse  and  the 
origin  of  the  coordinate  systems  are  at  the 
target's  estimated  position.   The  estimated  bear- 
ing a  and  estimated  range  p  are  indicated  for  a 
reference  position. 
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a2  =  — sy-i ~  Z  (1/a2)  (A  sin  3-  -  B  cos  3-)2 

J    (B2-AC)2  1      X  X  X 


and 


2. 


>_  =  — 5 £  (1/a .  )  (B  sin  3  .  -  C  cos  3  .  )  (A  sin  3  .  -  B  cos  3  •  )  . 

xy    (B  -AC)   1      x         1  1         x  1 


Using  the  definition  for  A,  B  and  C,  the  above  become 


tc\  2      c 

(6)  a.  =  


x    (AC-B2) 


(7)  a2  =  — ^ 


(AC-B2) 


and 


B 


(8)  ~~         2 

l0;  xy    (AC-B^) 


So  tan  2   =  2B/(A-C)  for  3-  =  <$>  ■  ,    i  =  1/  2,  ...,  n. 

With  the  target's  position  known  and,  consequently,  <j> .  known 

2    2 
for  i=l,  2,  ...,  n,  the  above  expressions  for  a-  ,  a~  ,a~~  and  v 

t-  x  '   y  *  xy      ' 

can  be  used,  since  the  initial  estimate  of  the  target's  position 

can  be  taken  as  the  target's  position. 

With  values  for  a~  ,   a-  ,  a~~  and  y,  values  for  o~,    and  o~, 

x    y    xy  x        y ' 

can  be  found  by  using  equations  (3)  and  (4).   The  probability 
that  an  estimated  position  will  be  within  an  ellipse  of  semiaxes 
kaA,  and  koA,  which  is  centered  on  the  target's  position  is 
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2 
1  -exp(-k  /2) .   This  result  can  be  found  by  integrating  the  bi- 

variate  normal  density  over  the  ellipse.   And  the  area  of  the 

2 
ellipse  is   k'  o~,o~,. 

x'  y 

Given  estimates  x  and  y  found  by  using  Equations  (1)  and  (2)  , 

the  ellipse  with  semi-axes  kaA,  and  ka~,  in  a  coordinate  system 

x         y 

that  is  centered  on  the  point  (x,y)  and  has  been  rotated  through 

2 

an  angle  y    is  a  1  -  exp(-k  /2)  confidence  region.   This  follows, 

since,  to  the  degree  of  the  approximations  involved,  the  bivar- 

iate  normal  distribution  of  X  and  Y  is  centered  on  the  target's 

.  .  2    2 

position.   The  confidence  ellipse  is  defined  if  a~  ,  a~  and  a~~ 

^  x    y       xy 

can  be  found,  that  is  if  the  elements  of  the  covariance  matrix 
can  be  found.   To  the  degree  of  the  approximations  involved, 
this  can  be  done  as  follows:   First,  assume  the  initial  estimate 

of  the  target's  position  is  at  the  target's  position.   Then, 

2    2  ... 

values  for  a~   a^   a~~  and  y  can  be  determined  by  using  Equations 
x  f  y  i    xy  J  ^ 

(6),  (7)  and  (8).   These  values  can  then  be  used  to  determine 

2      2 

a~,  ,  a-,  and  0*,~,    by  using  Equations  (3),  (4)  and  (5).   Now, 

with  a  value  for  k,  a  confidence  region  can  be  constructed.   To 
the  degree  of  the  approximations  involved,  the  shape  of  the  con- 
fidence region  is  independent  of  both  the  target's  position  and 
of  the  initial  estimate  of  the  target's  position. 

For  the  case  where  bearings  are  taken  from  the  target  on  two 
or  more  stations,  0.  is  the  reciprocal  of  the  bearing  taken  from 
the  target. 

A  discussion  of  the  theory  of  bearings  only  position 
estimation  procedures  for  situations  similar  to  the  one  considered 
here  is  given  in  Reference  1.   Reference  2  gives  an  essentially 
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equivalent  bearings  only  procedure.   It  also  gives  a  range  only 
procedure,  a  range  and  bearing  procedure  and  HP-9830A  programs 
with  which  to  implement  the  procedures.   Using  the  fix  determined 
by  two  lines  of  bearing  as  the  initial  estimate  was  suggested  by 
this  reference. 

The  following  equations  are  evaluated  in  the  program  to  deter- 
mine the  coordinates  of  the  initial  estimate: 

x*  sin  (92~9,)  =  [p-.  sin  (a, -6,)]  sin  02 


-  [p2  sin  (a?-62) ]  sin  6, 


and 


y*  sin  (92-9,)  =  [p,  sin  (a, -9,)]  cos  92 

-  [p2  sin  (a2-92) ]  cos  9,  . 

Reference  3  describes  a  TI-59  program  that  is  based  on  an  equi- 
valent procedure.   The  program  allows  a  user  to  either  input 
the  coordinates  of  the  initial  estimate  or  determine  them  in 
the  manner  described  here. 
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Appendix  2.   HP-41CV  Program  Labels 

The  global  labels  in  the  HP-41CV  program  that  are  assigned 
to  the  keys  E+,  1/x,  /x,  LOG  and  LN  give  these  keys  a  mnemonic 
character.   For  example,  with  the  calculator  in  USER  mode,  if 
LOG  is  pressed  and  held,  SIZ  will  be  displayed  and  then  after 
a  delay,  NULL. 

If  a  global  label  is  replaced  by  a  local  label,  this  mnemonic 
character  will  be  lost.   However,  automatic  key  assignment  will 
be  gained  if  the  key  corresponding  to  the  local  label  has  not 
been  previously  assigned.   If  this  is  the  case  and  if  there  is 
a  second  program  in  program  memory  that  uses  the  same  local  label, 
then  it  will  be  automatically  assigned  when  the  calculator  is 
positioned  in  program  memory  at  that  program. 

Global  labels  in  the  HP-41CV  program  can  be  replaced  with 
local  labels.   If  this  is  done  as  described  below,  the  user  in- 
structions will  still  be  applicable.   First,  make  the  following 
replacements: 

Line  329:  LBL  "PRB"  -*■  LBL  E.  Line  311:  LBL  "SIZ"  +  LBL  D. 
Line  186:  LBL  "EST"  ->  LBL  B.  Line  008:  LBL  "CON"  ->  LBL  C. 
Line  04  7:   GTO  "CON"  ->  GTO  C. 

Next,  after  line  001:   LBL  "TPE",  insert  LBL  A  so  that  it 
becomes  line  002.   Finally,  if  there  is  a  key  assignment  for  a 
key  in  the  following  list:   £  +  ,  I/x,  /x,  LOG  and  LN,  remove  it. 
Now,  with  the  calculator  in  USER  mode,  if  LOG  is  pressed  and  held, 
XEQ  D  will  be  displayed  and  then,  after  a  delay,  NULL.   The  key's 
mnemonic  character   has  been  lost,  but  if  there  is  a  second  pro- 
gram in  program  memory  with  the  local  label  D,  LOG  will  be  auto- 
matically assigned  to  D  when  the  calculator  is  positioned  in  pro- 
gram memory  at  that  program. 
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